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Abstract 

In 1963, Lieb and Liniger solved exactly a one dimensional model of bosons interacting by a 
repulsive 5-potential and calculated the ground state in the thermodynamic limit. In the present 
work, we extend this model to a potential of three 5-functions, one of them is repulsive and the 
other two are attractive, modeling some aspects of the interaction between atoms, and present an 
approximate solution for a dilute gas. In this limit, for low energy states, the results are found to 
be reduced to the ones of an effective Lieb Liniger model with an effective 5-function of strength 
c e ff and the regime of stability is identified. This may shed light on some aspects of interacting 
bosons. 
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I. INTRODUCTION 


The physics of Bose gases is a fascinating and complicated field of research. Since it 
involves a many-body problem, analytical results are rare and in some parameter regimes, 
one can use approximations to describe experimental systems with very good accuracy. For 
example, for a weakly interacting Bose gas, Mean-Field approximation can be used to reduce 
the many body Hamiltonian into a one-body non-linear Schrodinger Equation, the Gross- 
Pitaevskii Equation OT- In the opposite limit, a strongly interacting one dimensional Bose 
gas can be mapped into a gas of free Fermions (Tonks-Girardeau gas, see, for example, 
M). Exact solutions in other regimes are highly desired. 

Simple models like the Lieb-Liniger (LL) model, that may not have direct experimental 
realization, may alert us to unexpected physical phenomena that are overlooked when “rea¬ 
sonable approximations” are made and motivate experiments [6H8]. The model introduced 
in the present work is of this type. 

In their seminal work from 1963 [S], Lieb and Liniger managed to solve exactly a one 
dimensional model for interacting bosons. They considered the Schrodinger equation for N 
particles interacting via a 5-function potential 
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where Xj is the coordinate of the j-th particle and c is the amplitude of the 5 function. 
Making a Bethe ansatz [19] 
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where kp n are the k vectors obtained by the permutation P (where [P] is its parity) of the set 
ki ,..., kj\r. Lieb and Liniger wrote Bethe ansatz equations for the fc’s by imposing periodic 
boundary conditions on a ring of length L [TO , TTj, 
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These N coupled equations are solved numerically and the energy 
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is calculated for the ground state and the excitations mmm- 

In the present work, we study a simple model which takes into account the range of inter¬ 
particle interactions without giving up the mathematical simplicity. It is a generalization of 
the LL model [9] where in addition to the repulsion there is also attraction. It is defined by 
the Schrodinger equation for N interacting particles of mass m, 
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where the inter-particle interaction is modeled as a sum of three 5-functions: The central one 
is repulsive (co > 0) while the peripheral ones are attractive (q < 0). This model is inspired 
by the Van-der-Waals potential which has repulsive and attractive regimes. By adjusting 
the parameters Co,q,/ of (|5]), one can model scattering from many inter-particle potentials 


In section II, we present Bethe ansatz equations for two bosons interacting via three 5- 
functions interaction potential and in section III an approximation is introduced, that allows 
to extend the LL Bethe ansatz equations to an arbitrary number of particles. The ground 
state solution for the approximate equations is found in section IV. Section V specifies 
the parameters of the regime where the gas is stable. The results and their experimental 
relevance are discussed in section VI. 


II. BETHE ANSATZ EQUATIONS FOR TWO BOSONS INTERACTING VIA 
THREE 5-FUNCTIONS INTERACTION POTENTIAL 

We start by writing Bethe ansatz equations for a simple case where there are only two 
bosons. In this case, the equations are intuitive. 

Consider two bosons of mass m trapped on a ring of length L and interact according to 
(|5|. It is convenient to write the wave function ij) in terms of center of mass coordinate, 
rq = (xi + X 2 ) /2 and relative motion coordinate, r 2 = (aq — x 2 ) /2, 

(ri,r 2 ) = ~^=e lkiri (j) (r 2 ) (6) 

where k\ = 2ixn/L and n is an integer so that periodic boundary conditions are satisfied. At 
the center of mass frame of reference, k\ — 0 and the wavefunction of the relative motion, 
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</>(r* 2 ), satisfies the Schrodinger equation 

Yp cP 11 1 

-^-^2 + ^ c oH r 2 ) + ^ c iH r 2- l /2) + - Ci d(r 2 + l /2) <f) (r 2 ) = E<j> (r 2 ) (7) 

which can be written also as 
’ h 2 <9 2 

~2m ^r 2 + C °^ ^ + ° 1 ^ ^ r2 ~ ^ + Q< * ^ + ^ ^ ^ = 2Ei ^ ^ ' ( 8 ) 

As usual in such cases, the wave function takes a different functional form in each of the 
four intervals [— |], [—|,0], [0, |] and The result for 0 (r 2 ) is 
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where C is a normalization constant while Q 0 and Qi should be determined. They are easily 
determined for the 5-function interaction since the jump of the derivative at the locations 
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where r 2 = 0, c* = cq or r 2 = Z/2, c* = q. This results in two equations for Q 0 and Qi 
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k 2 should be determined to ensure periodic boundary conditions (/' (r 2 ) = (/' (r 2 + L/2). In 
addition, cf) (r 2 ) = cj> (—r 2 ) so that (/' (r 2 ) = —</' (—r 2 ). Therefore, in particular, </' (r 2 )| r ., =L / 4 
must vanish, leading to 
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Now, we return to coordinates x\,x 2 . For this purpose, we use the relations: k\ = 
ki + k 2 ^/2,k 2 = (ki -k^j/2,ri = (aq + x 2 ) /2 and r 2 = (xi — x 2 ) /2 resulting in 

hn + k 2 r 2 = k 1 x 1 + k 2 x 2 (15) 


and 


hn - k 2 r 2 = k 2 x i + kix 2 . 


(16) 


The function -0 of (J6]) takes the form 
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The periodic boundary condition 0 (aq + y,aq) = 0 (aq — %,x 2 ) results in 
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which are identical to (14) (under the assumption k\ = 0, namely, in the center of mass 
frame of reference). In the derivation we used the fact that aq — x 2 —> aq — x 2 + L involves 
rotation around the circle and consequently all the signs are changed. 


III. APPROXIMATE BETHE ANSATZ EQUATIONS FOR AN ARBITRARY 
NUMBER OF BOSONS 

The two particle solution cannot be simply generalized to an arbitrary number of particles 
since for small interparticle distances, 


|Xj — aq| < /, 


( 20 ) 


the sign function in equation corresponding to varies substantially. For small l, the effect 


of the regime (20) may be negligible as demonstrated in what follows. This is reasonable 
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for a dilute gas where l <C L/N. In such a situation, the LL solution is valid with the 
replacement ) — Q o + 2Qi cos ((kj — k s ) 1/2), leading to 
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The kj are distinct, namely, the wave function vanishes if kj = k s for s ^ j as was shown in 
the original work of LL [9]. 


In the region where inequalities (20) are not satisfied for any of the particle pairs, the sign 
functions in (|9]) are all equal. Therefore, in this regime, (21) is a solution with the spectrum 
(22). There is a Hamiltonian that is different from the original one, for which (21) and (22) 


are eigenfunctions and eigenvalues even if some of the inequalities (20) are satisfied. It is just 
defined by the eigenfunctions and eigenvalues. For 1 = 0, this Hamiltonian and the original 
one are identical. If the spectrum and the L 2 -norm of the eigenfunctions are continuous in 
l, the relative difference in the spectrum and the wavefunctions (in the L 2 -norm) goes to 
zero in the limit l —> 0. If they are also differentiable as a function of /, then the relative 
difference behaves as Nl/L. 

We show that for the low energy states, the 3-5 function system can be replaced by a 
system with one (5-function of strength c e ff- 

We assume 

(kj - k s ) l < 1 (23) 

for all wave vectors kj. In section V, we show that this limit is relevant for the ground state 
and low-lying excitations of a dilute gas since k max < (const) ^ is small. In the leading order 
in kjl, 
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the error is of the order Nl/L. Comparing (22) with (|3]), one finds that for small kjl, 
the behavior of the present problem is indeed similar to the one found for one (5-function 
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in the leading order in kjl. Eq. (25) is the main result of the present work, and it enables one 
to understand the physics of the three 5-functions interaction in terms of the one 5-function 
interaction. Of particular interest are situations where c e // is very different from Co + 2q (the 
total strength of interactions). In order to find such situations, we define the parameters 


r = ci /cq and x = mcol/h 2 and rewrite (25) as 
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This is a surprising result. It is instructive to analyze the behavior of c e /// (cq + 2c,), Eq. 


(26), as a function of x. We are interested in the regime x > 0 and —0.5 < r < 0. At x = 0, 


the derivative of (26) is negative and therefore the function decreases. At 

(1 + 2r) 


Xo = 


(29) 


it turns out that c e ff = 0 (even though c 0 + 2q ^ 0). Higher values of x result in negative 
values of c e ff, namely, the effective interaction is attractive (even though Co + 2 q > 0). 
Schematic description of c e ff/ (c 0 + 2q) is given in Fig. 1. 
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Figure 1: (Color online) Schematic description of c e ff/ (co + 2 q) as a function of x for r = ci/cq = 
— 0.25. The inset expands the region where c e ff changes its sign and the red dot is (xq,0). 

The result c e ff = 0 at x = Xq is verified numerically (see Fig. 2) and will be discussed in 
what follows. In the two particle case it is exact. For a related result see [m nsj. 


IV. GROUND STATE ENERGY 


In the previous section, we derived the approximate Bethe ansatz equations (22) for N 


bosons interacting by a three 5-function potential (|5]). The solution for these N coupled 
equations, (ki, /c 2 , • • •, fcjv), can be used to calculate the energy of the gas 
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In the ground state, \kj\ are minimal (but yet kj are different, as in the original work of LL 
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Lieb and Liniger [9j managed to calculate the ground state energy in the thermodynamic 
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We see that if l = 0, the ground state corresponds to the choice rij = j, (j = 1,... ,N). 
This is true also for l ^ 0, as long as 9 is a monotonic increasing function of k. To see this, 
assume kj > k rni then, by monotonicity of 0, 9 (. kj — k s ) >9 (k m — k s ) for all s, therefore 
G (kj) > G (k m ) and G (kj) is monotonic. Since 9 is an odd function, G (kj = 0) = 0. The 
kj for the ground state are the smallest possible in absolute value, hence, we choose rij = j 
for the ground state. Therefore, in the monotonic regime, 
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where 9' (k) = 89 (k)/dk and kj and k 3+ \ are adjacent wave numbers. Typically, 9 i 
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monotonic and (33) is justified at the regime where (23) holds (see Sec. V for more details). 
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Here, A is the Fermi momentum (this should not be confused with fermionic systems!). The 
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In these variables, p(z) = 9 (A z), the density of states is g (z) = p (A z), and Eqs. (35), (36) 
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How does one solve Eqs. (39) and (40)? First, it is necessary to choose values for «o, on and 
d. These values are related to the parameters of the Hamiltonian via the Fermi momen¬ 
tum A which is unknown at this stage. One should only keep in mind that ai/a^ = Ci/cq 
and therefore the ratio ai/ao does reflect the ratio between attraction and repulsion in the 


Hamiltonian. The integral equation (40) (with parameters a^,ai and d) can be solved nu¬ 


merically. The solution, g (z), should be substituted in (39) in order to find 70 . By repeating 
the above scheme for different parameters, it is possible to plot the dimensionless energy e as 
a function of the dimensionless interaction strengths q 0 and 7 ; and the dimensionless length 
d. for small values of d, the energy e depends only on the effective strength of interaction, 


that is (25) in dimensionless units, 


leff = 7o + 2 7* + 


aid (270 + 27 1 + cgd'y 0 + la 0 d^ 0 ) 
1 — baoaid 2 — aid 


(42) 


The significance of this effective strength of interaction is demonstrated in figure 2. In this 


figure, we present the solutions e(a 0 ,cq,<f) that were calculated by solving (39), (40) and 


(41). In Fig. 2(a), the energy is plotted as a function of the total interaction strength 70 + 27 ; 
and different choices of d are represented by different colors. It is clear that the effect of d is 
not negligible. Even at the regime d ^ 1, it is evident that the value of l has a strong effect 
on the ground state energy. Furthermore, even for a given value of d, the total interaction 
strength 70 + 27 ; (which is proportional to Co + 2 q) is not in one to one correspondence with 
the energy and therefore cannot be used to characterize the gas. Fig 2(b) shows that in the 


regime d 1, the energy indeed depends only on 7 e ff of (42). The results are consistent 


with (26). 
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(a) 


(b) 




Figure 2: (Color online) The dimensionless energy e of (41) as a function of dimensionless interaction 
strengths for q = — cq/4 (namely, r = —0.25) and 0 < «o < 30. (a) e as a function of 70 + 2 7 /. 


Different lines represent different choices of d of (38), from top to bottom: d = 0 (blue), d = 0.02 
(green), d = 0.04 (red), d = 0.06 (turquoise), d = 0.08 (purple). Points where the effective 
interaction is attractive were excluded from the figure (these were supposed to appear in the bottom 


purple curve in the regime x > xq = 2, see Eq. (29)), so that the highest value of x which does 
appear in the figure is x = 1.94 and it corresponds to the purple point (7.7,0.114). (b) The energy 


e of (a), plotted as a function of 7 e // (Eq. (42)). 


V. REGIME OF STABILITY AND DEFINITION OF DILUTE GAS 


The Bethe ansatz equations (22) and the effective interaction (25), are valid only where 


(23) is satisfied. Therefore, it is important to identify the regime where (kj — k s ) l -C 1. 
In the original LL model, the ground state energy and the values of Us are maximal for 


strong interactions, c —>■ 00 , where k n = ‘^fn and n' s are integers n = — y,..., y. Then, the 


N 


N 


maximal absolute value of k is k max = ^ and for all j, s, 


(kj — k s )l < 


2t tNI 


(43) 


For dilute gas, the inter-particle separation L/N is much larger then the interaction range l 


and (23) is satisfied. 


The same argument can be written for the three 5-functions interaction potential (J5]) . 
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If 9 of (32) is a monotonic increasing function of k, the ground state is given by rij = j. 


(j = 1, ■■■ ,N) and k 

max N • 

Let us analyze the function 9 (k) and identify the regime of parameters where it is 
monotonic, first, note that 9 (k) is monotonically increasing if and only if / (k) = 
j [Qo + 2 Qi cos (kl/2)] is monotonically increasing. For k —> 0, / (k) = —p^c e // and there¬ 
fore it is monotonically increasing as long as c e // > 0. Hence, if c e // > 0, there exist some 
k* (which depends on the parameters c 0 , q , Z and does not depend on L and TV since 6 is 
independent of these variables) such that for all k < k*, 6 (k) is monotonically increasing. 
For the ground state, the states with the smallest \kj\ are occupied, namely, rij = j with 
j = 1,..., TV, and 

G (. kj) < nN. (44) 

9 is an angle variable and therefore it is bounded (actually, for very small k, 9 = —i r). Hence 

N 

\kj\ < (const)— (45) 

Jb 

and can be made arbitrary small. Now, by increasing L (or decreasing TV), one may tune 
the value of k max such that the conditions 

k <r k* 

™max ^ rv 


and 


kmaxl 1 


are satisfied simultaneously, the regime (23) of dilute gas is achieved and our solution is 
correct up to a term of order Nl/L. 

For a dilute gas there is a range of parameters where c e ff > 0 and the solution is stable. 
There is also a range of parameters where c e // < 0 and the system is unstable. 


VI. SUMMARY AND DISCUSSION 

In this paper, we analyzed a one dimensional dilute Bose gas for an extension of the LL 
model defined by (|5]). By dilute gas, we mean that / -C L/N, that is, the effective size 
of a particle l (for example, the Van-der-Waals radius of an atom) is much smaller than 
the inter-particle distance. Using this assumption and the Bethe ansatz, we derived the 


approximate equations for the spectrum (18) ,(19), (22). In principle, these can be solved 
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numerically. For low energies in this situation \kjl\ -C 1 and the model can be approximated 


by a LL model with one 5-function of strength c e ff given by (25) and in dimensionless units 


by (42). The error of this approximation is of order Nl/L. This is a good approximation 
for the dilute gas. The effective interaction c e // depends on q and Co but also on the ratios 
between the characteristic potential energy scales, q/Z and co/Z, and the kinetic energy scale, 
h 2 /ml 2 , of a particle trapped in a well of length Z. 

Naively one would expect that for small kj, c e ff ~ Cq + 2q. It turns out to be correct for 
relatively weak interaction energy. For stronger interactions, c e ff becomes very small and 
even changes its sign (see Fig. 1). Note that this result holds also for the two particle case 
where it is exact. It is a surprising result, verified numerically in Fig. 2 and its experimental 
verification should be considered a challenge. The knowledge of c e // enables to calculate 
the ground state and the low excited states if the conditions for stability are satisfied. In 
section IV, the ground state is calculated in the thermodynamic limit for a dilute gas. In 
particular, it is demonstrated to depend on all parameters via c e ff. We have shown that 
for a dilute gas there is a regime of parameters where c e ff > 0 and therefore the system is 
stable. For other parameters, c e // < 0 and the dilute gas is unstable. In this regime, the 
results of H2H8] regarding dynamics of attractive gas might be realized. If the gas is not 
dilute, we cannot determine the stability of the system. This theoretical model enables to 
predict qualitative features of interacting bosons for realistic systems. 

The potential (J5]) can be realized, for example, in optical lattices [Hj with tight harmonic 
trapping along two perpendicular directions (E <C hco±) and almost flat potential along the 
third direction. The inter-particle interactions are in three dimensions and can be modeled 
by a “delta shell” potential 


V(r) = 


3cph 

4 rf n mui± 

cih 

2r'£ ut e ou tmuj_ 

0 


for r < r in 
for r out < r < r out + £ ou t 
otherwise 


(46) 


Where r in ,£ out —> 0. In a previous work [16j, we calculated the scattering length a and 
the effective range r e of the potential ( J46| ) (see App.A of [IB]), wrote a three dimensional 
Schrodinger equation and integrated it over two axes to obtain a one dimensional equation 
of the form (|5]). This leads to the relations 


Tout = \/Sl/2, 


(47) 
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a = 


(cq + 2 ci ) 


and 


r, = 


Ahw_ 

2cil 1 2 3 4 5 6 7 2 a 

a (cq + 2q) 3 


(48) 


(49) 


As seen from (25), for l = 0, c e ff is proportional to the scattering length. However, for 
l 7 ^ 0, c e ff cannot be expressed in terms of a and r e . Therefore, it motivates introducing an 
effective scattering length that dominates the spectrum. 

From an experimental point of view, it looks that c e // is the only quantity that one can 
measure in order to characterize the inter-particle interactions (because it determines the 
spectrum). Hence, it makes sense to define an effective scattering length 


a eff 


Ce f f 

Ahuj± 


(50) 


This scattering length, which includes corrections originating in the non-vanishing interac¬ 
tion range, is unique for one dimensional bosonic systems. 
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